A DISCRETIZED FOURIER ORTHOGONAL EXPANSION IN 
ORTHOGONAL POLYNOMIALS ON A CYLINDER 



JEREMY WADE 



Abstract. We study the convergence of a discretized Fourier orthogonal ex- 
pansion in orthogonal polynomials on B 2 X [—1, 1], where B 2 is the closed unit 
disk in R 2 . The discretized expansion uses a finite set of Radon projections 
and provides an algorithm for reconstructing three dimensional images in com- 
puted tomography. The Lebesgue constant is shown to be m (log(m + l)) 2 , 
and convergence is established for functions in C 2 (B 2 X [—1, 1]). 



for t e 1 and 9 £ S , the unit sphere in R d , is a classical topic with many 
applications to approximation theory and image reconstruction. In recent years, the 
role of Radon projections in topics dealing with orthogonal polynomials and Fourier 
orthogonal expansions in terms of orthogonal polynomials have been studied. These 
results rely on results put forward in two papers in by Marr [6] in 1974 and Logan 
and Shepp [7] in 1975, which uncovered a relationship between Radon projections 
and orthogonal polynomials on B 2 , the unit disk in R 2 . This relationship has 
led to many results in multi-dimensional approximation theory involving Radon 
projections. In 2005, Bojanov and Xu used this relationship to find a polynomial 
on B 2 that interpolates the Radon projections of a function taken on sets of parallel 
lines in directions given by equidistant angles along the unit circle [3], while in 2004, 
Bojanov studied the polynomial interpolating Radon projections on parallel lines 
in arbitrary directions [2]. In 1998, Petrushev extended the relationship in [7] 
and [B] to a connection between ridge polynomials on B d and higher-dimensional 
Radon projections [3]. In 2006, Xu introduced a discretized Fourier orthogonal 
expansion on B 2 , where the discrete data are finite Radon projections, and then 
studied the convergence properties of this expansion in the uniform norm [10j . 
This discretized expansion has application to image reconstruction in computerized 
tomography, where X-ray data corresponds to Radon projections. This discretized 
Fourier orthogonal expansion was generalized to the domain B d in |llj . 

In this paper, we study a discretization of the Fourier orthogonal expansion of 
functions on the cylinder B 2 x [—1,1] in terms of orthogonal polynomials with re- 
spect to the measure w(x,y,z) = (1 — z 2 ) -1 / 2 , where (x, y) € B 2 and z € [—1,1]. 



1. Introduction 



The Radon projection of function on M. d , defined by 




Date: June 12, 2009. 

2000 Mathematics Subject Classification. 41A10, 42B08, 41A63. 

Key words and phrases. Discrete expansions, Radon projections, orthogonal polynomials. 



2 



JEREMY WADE 



The discrete data are Radon transforms, taken on parallel disks, which are perpen- 
dicular to the axis of the cylinder. This discretized expansion is a linear operator 
on C(B 2 x [—1,1]), and the Lebesgue constant is shown to be m(log(m + l)) 2 . 
While the construction of this discretized expansion uses on the results from the 
discretized expansion on B 2 in [10] , the proof of the Lebesgue constant of the dis- 
cretized expansion utilizes a different approach. This is due to our definition of the 
degree of a polynomial being the total degree, which precludes us from expressing 
the orthogonal polynomial series on B 2 x [— 1 , 1] as a simple product of two or- 
thogonal polynomial scries on B 2 and [— 1 , 1] , and hence from obtaining a result 
trivially from previous results on these regions. In estimating the upper bound of 
the Lebesgue constant, generating functions for orthogonal polynomial series are 
used to separate the estimates on B 2 and [—1, 1]. This approach also provides an 
alternative proof to the original results on B 2 . The discretized algorithm has ap- 
plication in 3-d computerized tomography, as the discrete Radon projections can 
be interpreted as discrete X-ray data. 

This paper is organized as follows. In section 2, background material is presented, 
including some information on orthogonal polynomials and the Radon projection. 
The Fourier orthogonal expansion is briefly discussed in section 3, and the algorithm 
is presented in section 4. Since the derivation of the algorithm is also contained 
in [TU], sections 3 and 4 are brief. A detailed proof of the Lebesgue constant is 
contained in section 5. 



2.1. Radon projections. Let / be a function defined on R 2 . We define L(8,t), 
for 9 £ [0, 2tt] and t £ R, to be the line {(scos9 + t s'm9, ssin9 - t cos#) : sel}. 
The Radon projection of /, Hg(f;t) is defined to be 



By restricting the domain of / to B 2 , we only need to consider line segments 
in B 2 in the definition of 1Zg(f;t). For this reason, we define 1(6, t) to be the 
intersection of L(9;t) and B 2 ; that is, 

7(0; f) = {(s cos 9 + t sin 9,s sin 9-t cos 9) : \s\ <\]\-t 2 }, 

and re-define lZg(f; t) with the integral taken over 1(9; t) instead of L(9; t). 

For functions on B 2 x [—1,1], we will use Radon projections on parallel disks, 
perpendicular to the axis of the cylinder. For these functions, we define the Radon 
projection by 



where (x,y) £ B 2 and z £ [—1,1]. 

Radon's famous result states that a function can be completely reconstructed 
from its Radon projections, provided the function satisfies some modest decay con- 
ditions. Radon gave an explicit reconstruction formula, which required all the 
Radon projections of a function to be known. We instead seek to approximate a 
function using only a finite number of Radon projections. 



2. Preliminaries 




f(x, y) dxdy. 
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2.2. Orthogonal polynomials. Let ft be a subset of E rf with positive Lebesgue 
measure, and x = (xi,X2, ■ ■ ■ ,Xd) £ ^2. Let a — (ai,a2, ■ ■ ■ , ay) € Nq, and define 

x = x 1 x 2 ■ ■ ■ x d , |a| = ai + ot2 + ■ ■ ■ + ad- 

We define a polynomial P in d variables to be of degree n if it is of total degree n, 
which means P is of the form 

p(x) = c " x "> 

\a\<n 

where c a are real numbers, and at least one c a with |a| = n is not zero. 

Let uj be a positive measure on f2 with finite moments; that is, all integrals of 
monomials are finite. A polynomial P of degree n is orthogonal on Q with respect 
to uj if 

P{x)Q(x)duj{x) = 

whenever Q is a polynomial of degree less than n. 

We denote by V n (Sl;uj) the space of orthogonal polynomials of degree n on H 
with respect to the weight function w. In the case that w = 1, we simply write 
V n (Q). It is well known that dim(V„(fi)) = ( n+ d _1 ); see, for example, |t|. We 
define II„ to be space of polynomials of degree less than or equal to n. We may 
write II n = X)fe=o (^) > anc ^ ^ follows that dim II„ = (™ J rf ) . An orthogonal 
polynomial sequence is an ordered basis of orthogonal polynomials for the space of 
polynomials. 

On the domain £1 = [—1,1], the Chebyshev polynomials of the first and second 
kinds are denoted by 

™ / > „ / \ sin ((n + 1)9) 

T n (x) = cosn0, U n {x) = U . a ' \ 

smv 

respectively, with x — cos 9. The polynomials T n (x) are orthogonal with respect to 
the weight function uj (x) = ^(1— x 2 )' 1 ^ 2 , and the polynomials U n (x) are orthogonal 
with respect to the weight functions -(1 — x 2 ) 1 ^ 2 . We denote by T n the orthonormal 
Chebyshev polynomials of the first kind. They satisfy 

f (x) = T (x) , f n (x) = V2T n (x). 

The zeroes of T n (x) are 

(2.1) zi yU := cos(7;^ n ) := cos ( ^ l ^ n > i = 0,l,...n- 1, 

and the n-point Gaussian quadrature associated with the Chebyshev polynomials 
of the first kind is given by 

i r 1 , dx i ^ 



(2-2) -/ f(x)^==*-^f( Zl , n ). 

It is well known that equality is obtained in (|2.2[) if the degree of / is less than or 
equal to 2n — 1. 

The zeroes of U n (x) are 

(2.3) cos9j. n := cos — — ^, j = l,2,...n, 



4 



JEREMY WADE 



and the Gaussian quadrature associated with the Chebyshev polynomials of the 
second kind is given by 

(2.4) -/ f{x)y/l - a; 2 cfe « - V / (cos (0 i>n )) . 

Again, equality is achieved in (|2.4|) if / is a polynomial of degree 2n — 1 or less. 
For f2 = i? 2 , it was shown in [7] that the polynomials 

(2.5) U n (9 jjn ;x,y) :=U n (xco80j tn + ysm6j tn ), j = 0, l,...,n 

form an orthonormal basis for V„(i? 2 ,w) with respect to the measure uj(x,y) = — . 
This basis plays an important role in the construction of the discretized Fourier 
expansion on B 2 . The remarkable relationship between Radon projections and 
orthogonal polynomials, discovered by Marr in 1974, is given below. 

Theorem 2.1. Thm. 1] If P G V m (B 2 ), then 



n e (P;t) = 2 Vl f u m (t)P(cos8,sm9), 
m + 1 

where U m (t) is the Chebyshev polynomial of the second kind of degree m. 

On B 2 x [—1,1], define the weight function ujz(x,y,z) := ~(1 — z 2 ) -1 / 2 , with 
(x,y) £ B 2 and z € [—1, 1]. It follows from the above facts that, for (x,y) € B 2 
and z € [—1,1], the polynomials 

(2.6) Pj jk , n (x,y,z) :=U k (dj,k;x,y)T n -k(z) , fe = 0, 1, . . . n, j = 0, l,...fc 

form an orthonormal basis for V„(-B 2 x [—1, l];w^). This basis plays an important 
role in the construction of the discretized Fourier orthogonal expansion on B 2 and 
B 2 x [-1,1]. 



3. Fourier Orthogonal Expansion and Main Result 
Standard Hilbert space theory allows us to decompose L 2 (Vt;u>) as 

oo oo 

L 2 ((];.) = 0V fe (J];.), / = ]Tproj fc /, 

k=0 k=Q 

provided w is a positive measure with finite moments on J7. The Fourier partial 
sum of degree m, S n , is defined by 

n 

S n f = proj fc /. 

fc=0 

It is a standard result that S n f is the best approximation to / in II„ in the norm of 
L 2 (Q; u). We will use these standard results to construct our discretized expansion 
on the cylinder below. We first list the results on B 2 , and then move on to the 
cylinder. 
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3.1. Discretized Expansion on B 2 . The results in this sub-section are due to 
Xu and are proven in [10]. Recall the orthogonal polynomials U k (9j. k ; x, y) in (|2.5j) 
form an orthonormal basis for Vk{B 2 ). Using Marr's theorem [B], Xu was able to 
express the Fourier coefficients of the Fourier orthogonal expansion of / in this basis 
in terms of the Radon projections of /. 

Theorem 3.1. pH Prop. 3.1] Let m > and f € L 2 (B 2 ). For < k < 2m and 

< 3 < k, 

f{x,y)U k (9j,k;x,y)dxdy 

B 2 

= ^-tV- / n^(f;t),u k (t)dtu k (cos(9 Jtk -^)), 

where <j> v = ^ ff . 

As a result of this relationship, the orthogonal projection of a function / in 
L 2 (B 2 ) onto V k {B 2 ) can be written in terms of the Radon projections of /. 

Theorem 3.2. |1Q[ Thm. 3.2] For m > and k < 2m, the projection operator 
proj fc from L 2 {B 2 ) to Vk{B 2 ) can be written in the form 

. 2m -i 

projk/foy) = — — - V- / K^{f;t)U k {t)dt(k + l)U k ((t> u ;x,y). 
2m + 1 * — ' 7r / i 

As a consequence of this theorem, the Fourier partial sum of /, S^m/ '> can also 
be written in terms of the Radon projections of /, 

^ 2m 2m ^ „\ 

(3.1) S 2m f(x,y) = ^——J2J2- n^(f;t)U k (t)dt(k + l)U k (<t> u ;x,y). 
Recalling (|2.ip . the expression 

is a polynomial of degree 2m if / is a polynomial of degree 2m. Hence, by multi- 
plying and dividing by a factor of yl — t 2 in (|3.ip . and then using the 2m-point 
Gaussian quadrature rule given in (|2.4I) to replace the integral with a sum, a dis- 
cretized Fourier orthogonal expansion which preserves polynomials of degree less 
than or equal to 2m — 1 is obtained. This discretized expansion, Aami is given by 

^ 2m 2m 

-4am (/) 0, y) = (2 )2 E n <t>» ^ ; cos ^' 2 ™) 

* ' i/=0 3=1 



2 m 

X 
fc=0 

?2 



J^(fc + 1) sin((fc + l)dj,2m)U k ((t) u ;x, y). 



As an operator on C(B ), Aim has a Lebesgue constant of mlog(m + 1). The 
upper limit of the sum in v is chosen to be 2m to eliminate redundancy in the 
Radon data. With the choice of 2m, the Radon projections are taken along parallel 
lines in (2m + 1) directions, given by equally spaced points on the unit circle. If the 
Radon projections were taken along parallel lines in 2m directions, then for v < m, 



and <l> v - 2m 



2-KV 

TT + 4> n U = 4>u+m, 
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and the identity 

"^•7T+0„ (/; COs(02m+l-j,2m)) = K<t>„ (fi COs(0j )2m )) 

shows these two Radon projections are the same. In practical settings, more Radon 
projections are desirable, so we choose to take Radon projections in an odd number 
of directions. 

3.2. Discretized Expansion on B 2 x[— 1, 1]. The construction discretized Fourier 
orthogonal expansion uses the basis for V n {B 2 X [— 1, 1]; u>z) in (|2.6[) and the results 
on B 2 . 

Theorem 3.3. Let m > and let n < 2m. Define <p u := ^jj, an d a v{x,y) '■— 
arccos(x cos(</v) + ysm((f> v )). The operator proj„ can be written as 

(3.2) 

2m r i r i 



1 1 f f ds 

proj n f{x,y,z) = - 2m + l Ylj 1 J Kt v {f{-r,s)\t)y v {x,y,z\s,t)dt-j= 

1 " 

4V„0, y, z; s, t) = — VVfc + l)U k (t)U k (cos(a u (x, y)))f n _ k (s)f n _ k (z) 

2m + 1 



Proof. Since the polynomials Pj,k,n m (|2.6|) form an orthonormal basis for V n {B 2 

[-1,1]; w z ), 



proj„ /(at, y,z)=yyy\- 1 \ 

ZT n TTn 71 J-i Je 



n k 

f(x, y, z)Pj, k ,n{%, V, z)uz{x, V, z). 



After expanding Pj,k,n(x,y, z), Theorem fj3 . 2|) gives the result. □ 

This relationship between the projection operator and the Radon projection also 
yields a connection between the partial sum operator S2m and the Radon projection. 

Corollary 3.4. Let m > 0. The partial sum operator S^m may be written as 

2m 1 ^ * ds 



t) , dt 



(3.3) S 2m (f)(x,y,z) = V - / / n^(f(-,-,s);t)$ u (x,y,z;s, 
where 

2m 

§ v (x, y, z; s,t) = *^,„(a;, y, z; s, i). 
«=o 

We will discretize S^-m by approximating the integrals with Gaussian quadra- 
tures. In order to achieve a result on the convergence of our discretized expansion, 
we must first prove a result concerning the Radon projection of a polynomial on 
B a x[-l,l]. 

Lemma 3.5. If P is a polynomial of degree k on B 2 x [—1, 1] , then for 9 G [0, 2tt], 

n 9 {P{;;8)\t) 

VT^t 2 

is a polynomial of degree k in t. 
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Proof. If P is a polynomial of degree fc, we may write 

k 

P(x, y ) z)=y~" j CiZ % pk-i(x, y), 

where pk-i{x,y) is a polynomial of degree k — i in (x,y). Following the proof of 
[TOl Lem. 2.2], we write 

Ke(P(;;z);t) „ 1 



i=0 

Rewriting the integral and changing variables, 
p k -i{x,y)dxdy 

Wt) 

i-t' 2 



E Ci2;I ^=^ I Pk-i(x,y)dxdy. 
,-n VI - £ Ji(e:t) 



1 /-vx-t- 

= — / pfe_j(fcos^ + ssin#, tsin# — scos#) ds 

\/l - t 2 J-^T—p 

= J pk-i(tcosd + \J\ — t 2 s sin 9, t sin 6 — y/T— t 2 scos8) ds. 



After expanding pk-i (x, y) in the integrand, we note that each odd power of yl — t 2 
is accompanied by an odd power of s, which becomes after integrating. Hence, 



l r L ^ J 

== / Pk -i{x,y)dxdy = V b 3 t k -*- 2 Hl - t 2 y 



3=0 



for some coefficients bj. The lemma follows. □ 

With this lemma in mind, we substitue 

= — VI -t 2 

into p.3p . We then discretize the two integrals by using 2m-point Guassian quadra- 
tures. For the integral in s, the quadrature formula (|2.2p is used, and for the integral 
in t, we use the quadrature formula (|2 .4[) . The discretized Fourier expansion, S% m , 
is given below. 

Theorem 3.6. For m > 0, (x,y) € B 2 and z E [—1, 1], we define 

2m 2m 2m— 1 

(3.4) 6 2m (/)(x, y,z) :=EE E ^ W 

)) T v,j,z(a;,y,2:), 

v=0j=l 1=0 

where 

^ 2m n 

(3.5) T„jj(x, y, z) = E E( fc + X ) sin « fc + 

^ m ' n=0fc=0 

x J7 fc (cos(cr ly (x,j/)))T n _ fc (z ; )T n _ fe (z). 
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3.3. Main Theorems. As a result of lemma 13751 and the fact that 2m-pomt Gauss- 
ian quadratures are exact for polynomials of degrees up to Am — 1, we obtain the 
following theorem. 

Theorem 3.7. The operator Bi m preserves polynomials of degree less than or equal 
to 2m — 1; that is, for [x, y) G B 2 and z G [—1, 1] , 

B 2m (f)(x,y,z) = f(x,y,z) 

for f G n 2m _i. 

The proof of the following theorem is contained in the next section, and is the 
main substance of the paper. Comparing this with the result in [10], we see that 
the extension of the algorithm to the cylinder introduces a factor of log(m + 1). 

Theorem 3.8. For m > 0, the norm of the operator £>2 m on C(B 2 x [—1,1]) is 
given by 



||-B 2m ||oo ~ m (log(m + l)) 2 . 

The proof of this theorem is not trivial. Since we have defined the degree of a 
polynomial to be its total degree, the series in the definition of T Vt jj(x,y,z), 

2m n 

H( fc + X ) sin (( k + 1 )%.2m) U k (cob(ct„(x, y))) x f„_ fe (z / )f n _ fe (z), 

cannot be written as the product of a two series, one in terms of z and Zi , and one in 
terms of 9j and a v (x, y). As a result, the estimate of the Lebesgue constant cannot 
be trivially reduced to an estimate on B 2 and an estimate on [—1,1]. In particular, 
for the upper bound of the estimate, a different approach from that which was used 
in [10] is used to obtain our result. 

As a corollary of Theorems (|3.8p and (|3.7p , we obtain the following corollary. 

Corollary 3.9. For f G C 2 (B 2 x [—1,1]), 02m(/) converges to f in the uniform 
norm. 

Proof. If / e C 2 (B 2 x [—1,1]), then by Theorem 1 in pQ, there exists a polynomial 
p n of degree n on B 2 x [—1,1], and a constant C > 0, so that 

C ( \ 
||/-Pn||oo < -jWp - 
rr \ n 

where 



I \ 

sup \Wf(x) - W(y)\ 

B 2 x[-1, 



|J,J£B 2 X[-1,1] / 
\ lx-»|<l/n / 



|7|=2 

We let n = 2m — 1 to obtain 

||#2m(/) - /||oo < ||^2m(/ - P2m-l)||oo + ||/ P2rn— 1 1| c 

<ll/ 

P2rn— 1 1 1 00 (l + ||flam||oo) 
1 



< C- 



— (m(log(m + l)) 2 + l), 



(2m- 

which converges to zero as m approaches infinity. □ 
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Before proceeding to the proof of Theorem (|3.8|) . we make one comment. We 
believe that the Lebesgue constant of the Fourier partial sum of the orthogonal 
expansion, ||<S2m||j i s ~ mlog(m+l). We are able to prove HS^rolloo = 0(m\og(m + 
1)), but we have yet to prove the lower bound. If our intuition is correct, the 
discretization of the expansion adds a factor of log(m + 1). 



4. Proof of Theorem | 
We first derive an expression for which we may estimate ||B2ml|oo- 
Proposition 4.1. The norm of B^m as an operator on C(B 2 x [—1,1]) is given by 

2m 2m 2m— 1 

(4.1) ||S 2m ||oo = 2max^ ^ ^ sin0j, 2m \T u<jt i(x, y, z)\ 

v=0j=l 1=0 

where the maximum is taken over all points (x,y,z) in B 2 x [—1, 1]. 
Proof. By definition, 

(4.2) (/(-,-, zz), cos 0j,2m) = / f(x,y,zi)dxdy 

J / (COS #j,2tti ,<^v) 

sin Oj^m 

/(cos 9j.2m cos 4> v — s cos (\) v , cos 9j.2 m sin 4> v + s cos 4> v , z{) ds. 

- sin 9j i2 m 

Taking absolute value of both sides and using the triangle inequality, we immedi- 
ately have 

Ira Ira 2m— 1 

ll^mlloo < 2max^ ^2 ^2 sin9j,2m\T U j,i(x,y,z)\ 

v=0 j=l 1=0 

On the other hand, if we define 

Ira Ira 1ra—\ 

T(x,y,z) :=2^^ sm0j,2m \T„,j,i(x,y,z)\ , 

v=0j=l 1=0 

then T(x, y, z) is a continuous function on B 2 x [—1, 1], and hence achieves its maxi- 
mum at some point (xo, yo, zo) on the cylinder. We would like to choose a function / 
so that f(x,y,z) = sign(T v j t i(xo,yo,zo)) on the set of lines {(/(cos#j,2m, <t> v ), zi)}, 
for 1 < j < 2m, < v < 2m, and < I < 2m — 1, since this would immediately 
give us the result. However, such a function may not be continuous at the points of 
intersection of these lines. To allow for continuity, we instead take neighborhoods 
of volume e around each point of intersection of the lines, and define a function /* 
which is equal to sign(T yj -,i{xo, yo, z )) on the lines {(/(cos 6^2™, <f> v ), except 
on the e-neighborhoods at the points of intersection; on the rest of the cylinder, 
/* is chosen so that it takes values between 1 and —1 and is continuous. It then 
follows that 

Ira 2m 2m— 1 

ll^mlL > \B2m{f*(xo,yo,Zo))\ > 2 X! Sm ^,2m (x , ?/0 , ^o) | ~ CE, 

v=0j=l 1=0 

where c denotes the number of points of intersection of the lines {(I(cos9j.2m, 4>v), Zl)}j,v,l- 
Since e is arbitrary, this proves the proposition. □ 
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For the remainder of the proof, the number n in (|2.1|) and (|2.3| will be fixed as 
2m. For this reason, we define 



(4.3) 



0,2m _ 

2irv 



JTT 



2m +1 



21 + 1 
Am 



-IT, Zl = Z lt2m = COS 7/, 



o u (x,y) — arccos(x cos 0„ + ysm.<p v ). 



" 2m+f 

The proof will be separated into two parts, a lower and an upper bound. 



4.1. Lower Bound. We will establish there exists a constant c > so that 
||£>2m||oo > cm(\og(m + l)) 2 for all m > 0. By (|4.ip . it suffices to show 



2m 2m 2m— 1 



X! sm ^' I^.K^yi^i)! > cim(log(m + 1))" 

i/=0j=l 1=0 



point (xi, j/i, 2i) = (cos- 



, sin ■ 



, 1 ) for some c > 0. We begin 



for the jj^ii^u \u, L ,y U A, L j — I ^ 4m+2'"" 4m+2' 

by deriving a compact formula for T u j t i(xi,yi } 1). Using the Christoffel-Darboux 
formula for T„, letting cos7 z = z, 



\Tv,j,i(x,y,z)\ = 



1 



(2m + l) s 



2 m 



V sin((fc + l)ov(x, y)) sin((fc + 1)0 j) 

sm fj^^x, yj 



k=0 



cos((2m — fe + 1)7 Z ) cos((2m — k)ji) — cos((2m — k + cos((2m — k)j z ) 



cos(7 z ) - cos(7 ; ). 

Substituting in z = 1 and applying the identity for the difference of cosines, 
\T„,j,i(x,y, 1)| 



(2m + l) 3 sinov(x, y) sin : 



2 m 



^2(k + 1) sin((fc + l)a u (x, y)) sin((fc + 1)0,-) sin ((2m - fc + 1/2)7,) 



fe=0 



Applying the product formula for sine and the product formula for sine and cosine, 



(4.4) |T^(x,y,l)| = - 



1 1 



2m 



k=0 



4 (2m + l) 3 sin ^ sma u (x, y) 
J2(k + 1) [sin((fc + l)(0j - a v {x, y) + 7i ) - 

- sin((fc + l)(9j - a v {x,y) - 7z ) + | 7; + f ) 

- sin((fc + l)(0j + a v {x, y) + 7i) - §7; - f ) 

+ sin((fc + l)(0 3 +<r 1/ (x,y) -7i) + §7; + f)] 
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Next, apply the formula 

(k + 1) ain((A! + 1)0 + 4>) 

1 (2m + 2) sin((2r7i + 1)6* + 4>) - (2m + 1) sin((2m + 2)6 + (f>) + sin(</>) 



2m 

E 

fc=0 



1 sin((2m + 1)6 
2 



Bin 2 (I) 

(4m + 2) cos((2m + 3/2)0 + 0) sin(0/2) + sin(^) 



sin 2 (f ) 

to (|4.4p . Under our choice of X\ and cos av (x i,yi) — cos 2 ^r^ 7T) so cr v {x\,yi) = 
2 2m+i 77 ^ ^ — ^ — m ' ^ e WU ^ on ly De considering v within this range, so we define 

2i/-l/2 „ 
2m+l 



(-ly+'+i C o S (a) + (_l)J+'+l(4m + 2) sin ($±*) sin (tt^j ± cos ) 



x 2 / 9+0+7 



Taking into account 

(2m + l)(0i,2m ±1 ov ±2 70 = (i ±1 (2^ - 1/2) ±2 (I + 1/2))tt ±2 7/, 

where the subscripts indicate the signs of ±1 and ±2 are not related (a convention 
we will adopt for the remainder of the paper), we are able to write 

(4.5) 



',5,1 ( cos 



4m+2 



, sin 



4m+2 



,1 



1 



(2m + l) 3 8sin<r„ sin y 



fZitfj, H) ~ F^tTu, li) - l,j!<>,- Tr. ~li) + F+(6 j)( T u , - 7 i) 



T, 



We will show the lower bound is attained if we restrict the summation in (|4JJ) to 
the set of indices where 7r/4 < 9j < 3ir/8, ir/4 < 7; < 9j, and < er„ < 9j — 7/, so 
we only take the sums over the following range of indices: 

• LfJ+5<i<3LfJ 

• Li J + 1 < 1 < j - 4 

• 1 < v < - 1 

We assume m > 24, so that these inequalities make sense. With this restriction of 
summation, sin0j and sin y are bounded away from zero by a positive constant. 
Hence, we are left with proving the estimate 

(4.6) 

3LxJ i-4 L^J-i 

E E E ■ 

i7 = L^J+5i=L^J+l *=i 



1 



1 



(2m + l) 3 sinov 

x J F t & ' Tl) (0, , o v , <yj) - F u (9 3 ,-* u ,- 7l ) + F+ (9 3 ,a v ,- 7 i) 
> cm(log(m + l)) 2 
Also note that, under this restriction of summation, 



o< 



7i 3 

< -7T, 
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SO 

; ±1 oy ± 2 "fi 
2 

where we have used the fact 



±1 ov ±2 lu 



sin 9 w 6» 

if — 157r/16 < 9 < 157r/16, a fact we will use repeatedly throughout the proof. 
The dominating terms in the summation will be the terms 



(4m + 2) an 
sin I - 3 



2 

lemmas to eliminate the non-dominating terms. The first lemma eliminates the 

3,1 

Lemma 4.2. Recalling (|4. 3[) . 



the middle term in the numerator of from i(9j, ±cr„, — 7;). We first prove two 

rm 

first and third terms in the numerators of f - j 1 (0j, ±20V, ±37;). 



7=^* E E E 



(2m+l) 3 ^ ^ ^ a ir,rr a ^r, 2 f 2i±l2j^22l\ 

v ^ ; J= Lf j+5/=Lf J+i smff " sm ^ 2 ) 

< cmlog(m). 

Proof. First, considering 9j ± er„ + 7;, apply the inequalities 

% + ov + 7; > Oj + 7/, and 9j - a v + 7/ > 3 - + 7/ - 7r/8, 

to obtain 

x 3 LtJ i-4 j L^J-i x 

Jl(±1 ' +) - (2m + 1) 3 . ^ 5 cin 2 f te^ T ^ srn^: 

- c E p E E 21/ - 1/2 

^Lfj+s- 7 i=LfJ+i "=1 7 

< cmlog(m + 1). 

For J(+, — ), using the inequality 9j + a v — 7; > 0j — 7;, 

j 3Lfj i_ 4 x L*rJ-i x 



(2m + 1) 3 4-; , 2 / 8,-7^ ^ sin(ov) 

^ c E E ,.. ; ^2 E 



(i-l-2) 2 ^ 2^ — 1/2 
< cm log(m + 1). 
For the remaining case of J{—, — ), we split the sum in v, 



E + E 

f=i v=yi^L\) su^oy) sin 2 



1 



2 / Bi-o v -n 
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We are only considering values of v > 1, so we ignore any instances of v = 0, — 1 in 
the sums. For the first sum, 6j — er„ — 7, > (0j — 7/)/2, so a bound of cm log(m + 1) 
is found as in the case of J(+, — ). For the second sum, 

21/ -1/2 > (j-l- 3)/2 > 1 
2m + 1 " 



so it readily follows that 

3L^J j-4, 

*(-,-)<« E E 

j=Lf J+5(=Lf J+i 
+ cm log(m + 1) 
< cm log(m + 1). 



2m + 1 ~ 4m + 2 ' 



E 



J- l ~3 t~[ (j-2i/-f-l 



4m / 



□ 



The next lemma eliminates the parts of F- z (0, ^fa v , 7/) with 4m + 2 in the nu- 



merator. 

Lemma 4.3. Recalling (|4.3[) , 



■M±) := 



j_ 3L^j ,-4 ^-1(401 + 2) sin 

3 E E E 



(2m + 1) 



3 =[f J+5(=Lf J+i "=i sin ov sin 



< cmlog(m + 1). 



Proof. Since 3tt/4 > 0J+OV+7, > d. } -G v > tt/4, both sin (^p 1 ) and sin ( gj± ^ +7 ' ) 
are bounded away from zero by a positive constant. The lemma then follows from 



3L^J j-i L^J-i . 

— E E E 5^ 

i=L^J+5J=L^J+i f=i 



2m 

< cmlog(m + 1). 

By applying the triangle inequality to (|4.6[) . we obtain 

1 i=i 



-1/2 



□ 



IB; 



oc 



> 



1 



(2m + l) 3 



3L^J LVJ-i , 

V V V — 



j=L^J+i/=Lf J+i "=1 v y 
(1^(^,^,-71) -^7,(^,-^,-71)1 - 1^,(^,-^,71)1 - 1^7,(^,^,71)1) • 



The two lemmas show that 

3L^J j-2 L^J-i 



1 



(2m + l) 3 



E E E ^(^(^-^^l + l^^^^l) 



< cm log(m + 1). 
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We also have 



FU 6 ^ -V) - F u( 6 ^-^ "TO > (4m + 2) 



sin 



(M^) sin(^) 



sm 



sm 



Lemma (14.21) shows that 



2 f gj-QV-Tj ^ „;„2 / e 3 +<r„-7i 



1 3L^j j-2 L^J-i j 

< cm log(m + 1). 
Hence, we obtain 
(4.7) 



\\B2mW00 > 



1 



(2m + l) 2 



3LtJ 3-2 L^J-i 

E E E 



J= L^j+5/=Lf J+i " 
cmlog(m +1). 



—J sin (ct„) 



,in(to) sin(^) 



We now show the dominant part achieves the bound of cm(log(m + l)) 2 . Using 
the formula for the product of sines, 



sin(^) sin (*±^) 



sin(te-) sin(to±-) 

and it follows that 

3LxJ 3-4 L^J-i 

E E E 



sn^ov) sin (^-) 



sin(te-)sin(to±-) 



1 t^i sin ( ff -) 



sm 



j=Lf J+i'=LfJ + 

3LfJ j-i L^J-i ! 

~ C E EE Tq _ 7 -a v ){6i -7; + oyj 
i=Lf J+i;=Lf J+i "=1 VJ AJ ' 7 

- E E ^ E 

J= Lf J+i/=Lf J+i 

3LfJ 3-4 

> c E E 



?i - 7/ 

L^J-i 



1 



1 



1 



E 



j- - 7; - er„ &j - 7/ + er„ 
1 
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Dividing by (2m + 1) , we have 

E E j^- E 



(2m + l) 2 ^ ^ 0,- - 7; ^ 0,- - ov - 7; 





i-4 


^ E 


E 






3LtJ 




^ E 


E 


i=Lf J+ii= 


= Lf J+i 


3LxJ 


J-4 


^ E 


E 




=L?J+i 


> cm(log(m 4 


-I)) 2 , 



3 =Lf J+i'=Lf J 

I V 



L^J-i 

E • 



■2ro+l ' 



log(j-0 
i - I - 1/2 



which completes the proof of the lower bound. 

4.2. Upper Bound. As mentioned previously, the estimate of the Lebesgue con- 
stant on B 2 x [-1, 1] does not trivially reduce to estimates on B 2 and [—1,1]. More- 
over, a straight-forward estimate of the upper bound of the Lebesgue constant, as 
done in |10j for the case on _B 2 , would be extremely lengthy and cumbersome. We 
instead use a different approach, by deriving generating functions for the series in 
the definition of T u j t i(x,y, z), and then writing T u j t i(x,y, z) as the Fourier coeffi- 
cient of the product of the generating functions. The idea for this approach comes 
from [5], and provides an alternative proof for the upper bound of [10l Theorem 
5.2]. 

Lemma 4.4. For < \r\ < 1 and m > 0, 

T v ,j,i{x,y,z) 



(2m + l) 3 siner„(:zi, y) 

T- 1^ r ^G 1 (re^,z,z i )G 2 (re l9 ,^,a 1 ,(x, y ))e- 2 "" e ^r- 2 '», 
2n J q 1 — re w 

where 

oo 

Gi(r,js,js,) :=J2 f k{z)f k {zi)r k , 

k=Q 

and 

oo 

G 2 (r, 6j,a v (x, y)) := ^> + 1) s H(k + l)(9j)) sm((fc + l)(a v (x, y)))r k . 
Proof. Our first step is to derive the generating function of the function 



RN{0j,cr u (x,y),z,zi) 

N n 

<j v {x, y))T n - k (zi)T n - k (z) 



ji=0 k=0 
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in T„ t jj,(x, y, z). Since the coefficient of r in 

oo n oo 

E E sin « fc + sin « fc + l ) a »( x > y)> n E f Az)TMy 

n=0 fe=0 j=0 

is precisely Rn(0j, a u (x, y), z, z{), and 
1 

— E + 1)^ ) sin((fc + l)a v (x, y))r k 

oo n 

= E E sin (( fc + !)^) sin (( fc + y)y 



r 

k=0 



n=0 fc=0 

it follows that 



E RN{9j,Vi>{x,y)re'z,zi)r k = ^—;Gi{r,z,zi)G 2 {r,6j,o- v (x,y)). 



N=0 



Since both sides of the above equation are analytic functions of r for |r| < 1, we 
may replace r with re* 6 * to obtain analytic, complex-valued functions of r. Since 



R2m(e j ,a„(x,y),z,zi) = — Y] R N {6 3 ,o- v {x,y) 1 z 1 z l )r 2m e 2mie e~ 2 ""■'»■■ 2 '" 

2n J o ^0 



it follows that 

R2m{0j,(J v (x,y),Z,Zl) 

= T" / 1 ^G 1 (re ie ,z,z l )G 2 (re ie ,e j ,a v (x,y))e- 2mie de r - 2m . 

The lemma follows from the fact that 

(2m + l) 3 sinc^(x,y)i? 2m (6^,av(2;,7/),z,z ; ) = T„jj(x, y, z). 

□ 

We next obtain compact formulas for G\ (re t6 , z, z\) and G^fVe* 6 ', Oj, <r u (x, y)), 
and obtain estimates for these functions. 

Lemma 4.5. For m > 0, and r = 1 — -, 

— ; m ' 

1 1 1 f 2 " 1 

(4.8) |T^(x,y,z)|< (2m + i) 3 s . n ^ (x ^ ) -y o 

x ( | (x, y) - (jc, y) | + | A+ (x, y) - (x, y) | ) 
x (\P(re w , lz + 7i ) | + |P(re^ 7z - 7 ,)|) d$. 

where 

(4 ' 9) = l-2re^cos{6 J ± ( j u {x, y )) + r 2 e^ 1 

and 

(4.10) ^ - ^) { re-- cos (*, ± „ „(*,„))) 



(1 - 2re ie cos (^ ± o v {x, y)) + r 2 e 2W ) 2 ' 
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Proof. First, it follows from Lemma (|4.4j) that 
1 1 



\T U!j j(x,y,z)\ < 



(2m + l) 3 sin a u (x, y) 



X ^ I h 1 . t9 , |Gi(re l9 ; z ; zO| |G2(r e l9 ,^-,^(x^))|^r- 2 '». 
2tt Jo I 1 - re I 

We next derive estimates for |Gi(re ie , z, z;)| and \G2{re i6 9j, a v {x, y)\. The compact 
formula for the generating function G\{re ie , z, zi) is well-known, 

(4.11) Gi(re» *,*,) = J [P(re* 9 , 7z + 7/) + P(re w , l z - 7l)] > 
where P{r, (f>) is the Poisson kernel, defined by 

00 1 - r 2 

(4.12) P(r,0) :=l + 2]Tcos(n0)r n = - — -, < r < 1. 

' 1 — 2r cos <p + H 

n— 1 r 

For G 2 (re , 9j,a y {x, y)), we use the identity for the product of sines to obtain 
G 2 (re ie ,^,0 
1 (7 00 

= Ydr £ t cos (( fc )^ - ^(^ »))) ~ cos (( fc + !)(^ + °»&V))\ ***** 

k=l 

= ~ [P(re i9 A- -<7„(!c,y)) -P(re ie ,^ + <7„(a;,y))] . 

Using the formula 

-P(r,0 = -2 f - | (^-co^Kl-r 2 ) 

dr ' \1 — 2rcos0 + r 2 (1 — 2r cos + r 2 ) 2 

we obtain 

re ie [re w - cos(0j - a v {x,y))](l - r 2 e 2ie ) 



1 - 2re lt > cos(0j - ov(a:, y)) + r 2 e 2l ° [\ _ 2re l6 coa(9 3 - a„(x, y)) + r 2 e™f 
re ie [re 19 - cos(6» j + a v {x,y))]{\ - r 2 e 2ie ) 



1 - 2re ie cos(9j + a v (x, y)) + r 2 e 2W [l - 2re ie cos(6» J + <r u (x, y)) + r 2 e™\ 
It follows that 

I G 2 (re lS , 6 j , a v (x, y)) \ < ( | A+ (x, y) - (x, y) \ + \ A+ (x, y) - (x, y) \ ) . 



Finally, since the inequality in (J4T4J) holds for all values of r with < r < 1, we set 
r = 1 — - . With this choice of r, r~ 2m converges to e~ 2 as m approaches infinity, 
and so r~ 2m is bounded by a constant for all m. □ 

Before beginning the estimate, we make several reductions in the range of the 
sums and values of x, y, z that need to be considered. 

(1) First, we can reduce the interval of integration to [0, ir]. To see this, replace 
9 with 2ir — 6. This change of variable amounts to conjugation of the 
complex number re l ° , and hence the norms of the expression are unchanged. 
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(2) We may also restrict j z to the interval [0,7r/2m]. To see this, replace 7 2 
with 7z + in P{re l9 ,-y z + 9j) and P{re lB ,-y z — 7/). We see that, upon 
changing the summation index from / to 2m — 1 — I, 

2m-l 



1=0 

2m-l 



v 1 - r 2 e 2lS 

f-f 1 - 2re w cos (j z + ^7 + 7) + r 2 " 2 '" 



1=1 
and 

2m-l 



2re ie cos ( 7z + 7; ) + r 2 e 2i0 + 1 - 2re l6 cos ( 7z - ^-n) + r 2 e 2i6 ' 



1 — r 2 e 

^ 1 - 2re if> cos ( 7z + ^ - 7) + r 2 e 2W 



2 Z^Z 2 l _ r 2 e 2 t 9 ! _ r 2 e2l e 



2re if > cos - 7) + r 2 e 2ie 1 - 2re i9 cos (j z + ^) + r 2 e 2if? ' 



It follows that the expression 

2m-l 

^ |P(re l9 ,7 z +7)+P(re ie ,7 z -7)| 

1=0 

is invariant under translations of 7 Z by 7r/2m, so we only need to consider 

(3) The sum in j may be reduced to 1 < j < m. Replacing j with 2m + 1 — j, 
sin^m+i-j = sin0,-, and cos(0 2m .+i-j ±v v {%, y)) = cos(6j=f(ir - a u (x,y))). 
It follows from the definition of a v {x,y) that 7r — a u (x,y) = a v {—x,—y), 
which implies that sin er„(— x, — y) = $mo v (x,y). Hence, 

2m s' 

y, sm (\A+(x,y)-Al(x,y)\ + \A+(x, y) - A^{x, y)\) 

j=m+l uy ,y ' 

m . ^ 

sm Uj 

sinov(-x,-y) 

x (-«, -V) - A l(-x, ~y)\ + \At(-x, -y) - A^i-x, -y)\) , 

which shows we only need to consider 1 < j < m. 

(4) We also only need to consider (a:, y) in the region 

(pcoB^prin.,) : "4^<^<^: 

To see this, let x — p cos 77 and y = psin#, so that cos ov(x, y) — pcos(rj — 
<f) v ). Note that the collection of points pcos(r] — 4>u), for u = 0, 1, . . . ,2m, is 
unchanged by a rotation of 77 by 0„. Moreover, every expression involving 
<j v {x,y) in the right side of (14. 8| can be written in terms of cosov(:z:, y). 
Since < a v {x,y) < tt, sma u (x,y) = y/l — cos 2 av(x, y), and the expres- 
sions cos(6j ± ov(a;, y)) can be expanded using the cosine addition identity. 
Hence, every expression involving a u (x, y) in 

2m 

Y,-. 7 r - Ai (as.y)! + |A+(x, y) - A^(x, y)|) 
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is the same at the points (p cos 77, p sin 77) and (pcos(j] + 0„), psin(?y + <f> v ))- 
(5) Finally, we may also reduce the sum in v to < v < m. First note that 

cosa 2m +i-v(x,y) = cos a v (x, -y), and sina 2m +i- v (x, y) = sinov(:r, -y). 

Hence, we obtain 

2m 

E sinfT ( x (l^^'f) - A i( x >y)\ + \ A i( x >y) - A 2(x,y)\) 

y—ra-{-l V ^ ' ' 

m _^ 

= E - — ? r (l A i~( x > -f) - + 1^2 - ^fo , 

^j— J srno-„(a;, -y) 

and since r m is symmetric with respect to y, we only need to consider 
< v < 771. 

From these reductions, it follows that 

r 2ir -t 2m— 1 



|S 2m ||oo < c max / — V \P( 

(x.,y)&T m J Q 2m + 1 ^ 



^ y ,7 2 +7i) + ^(^,7,-7i) 



.| 2 m 2m 

EE' 



(2m + l) 2 |1 — re t6 \ sma u (x, y) 

x (\A+(x,y)-Ai(x,y)\ + \A+ (x,y) - A^ (x,y)\) d9. 
Consequently, the proof of Theorem l3 . 81 will follow from the three lemmas below. 



□ 



Lemma 4.6. For z £ [0,ir/2m], 

2m- 1 

2m + 



^ 2m- 1 

— |-P(^e ie ,7 2 +7/)| + |^e ie ,7 z -7/)| <clog(m+l) 



1=0 

for some c which is independent of 9, z, and m. 
Lemma 4.7. For (x,y) £ T m , 

-t pit 1 2m 2m . ,„ \ 

/or some c which is independent of x,y, and m. 
Lemma 4.8. For (x,y) € T m , 



TI f E E Si ! l(g f m l tf) - AZ(z, y)\d9< cm(log(m + 1)) 



(2m+l) 2 y |l-re ie | sin^Or, y)) 

for some c which is independent of x,y, and m. 

The proofs of these lemmas are contained in the next section. 
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5. Proofs of Lemmas [4761 14.71 and 14.81 
5.1. Proof of Lemma 14.61 

Proof. Fix z € [0, . First applying the triangle inequality, we have 

|Gi(r,MI 
< c 



1 — r z e 



2„2j<9 



1 — r z e 



2„2i8 



1 — 2re l9 cos(7 z — Ji) + r 2 e 2ie 



|l - re l0 \ X | sin § | + m" 1 



1 — 2re lS cos(7 2 +7;) + r 2 e 
Next we apply the approximation 
(5.1) 

and also note that 

(5.2) 1 - 2re lB cos(0) + r 2 e 2l9] = (1 - re l(e+0) )(l - re 1 ^"^) 

to obtain 

|Gi(r,0,z)l <c 



(| s in(^ 


) +TO" 1 ) 


;|sin( 


§)|+m- 




( 


sin( e+ ^+^) 


+ m-i) ( 


sin(* 


2 J 


+ m _1 ) 





(|sin(^ 


) + m" 1 ) 


J sin (|) + m 




( 


sin ^ 


9+72-71 ^ 
2 ) 




S in(^) 





:= c(*i + * 2 ) 

Case 1: < 9 < tt/2. 

We may ignore the factor of sin f 2 ^-) in this case. We also have 



1 + 9 Z + li 



< 



3~K 7T 

4 4m' 



so that sin((# ±i j z ± 2 7z)/2) ±i 7 2 ±2 ji- Fixing 9, for vE^i, we have 

2m-l 2m-l 



1 



1 ^ 2m + 1 ^ 

1=0 

2m- 



9 + rn 



2m + 1 ^ ^"2m + l ^ (0 + 7z + 7i + m" 1 )^ - 7z - 7z | + m" 1 ) 

2m-l 



2m + 1 f-< Q 10-7, -7i|+m- 1 

2m 



2m ^ 

< c > < clog(m + 1). 

(=0 

This type of estimate will be very common throughout the proof, and we will 
frequently omit the repetitive details. 
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For ^>2, we first sum starting from I = 1, so that 6 + "/i ~ j z > 0, and hence 

i ^ 



2m +1 f-^ {\9 + lz ~ 1 i\+m- 1 )((\9- lz + ll \+ m - 1 ) 



2m- 1 

I 

< 



- — 7 ^ 77 — i r < clogfm + f). 



2m + 1 ^ |# + 7.z - 7;| + to 

For the term corresponding to / = 0, either -f z > 7t/4to or 7 Z < 7r/4m; assuming, 
without loss of generality, the former, it then follows that 

1 6 + m- 1 

2 m + i(\e + lz -^\+ m ^)(\e- lz + ^\+ m -^) 
< * L_ <!, 

- 2to + 1 \0- 7z + jl.\ + m -i " 
which shows that 

^ 2m- 1 

$> 2 < C l0g(TO+l) 
Z=0 

for the case when 9 € [0, n/2]. 

Case 2: 7r/2 < 9 < n. For this case, the factor of sin(0/2) can be ignored. For 
note that 

TT 2n - 9 - J z - Jl < 37T 7T 7T < 6> - 7 Z - 7; ^ 7T 

~ 4m~ - 2 - T' ~4 ~ 8m" ~ 2 ~ 2 ' 

so we may approximate the sine functions accordingly and obtain 

2m— 1 2m— 1 , «, , _i 

1 , c x - \n-9\+m 

1 ^ 2to + 1 f-^ 



2to+1 ^ "2to+1 ^ (\2n - 9 - lz - lt \ + m- 1 ) (\9 - lz - -fi\ + m- 1 ) 

2m— 1 I /II i — 1 

C x - |7T — 9\ + TO 

~ 2to+ 1 ^ (|tt - 6» - lz + 7; | + to- 1 ) (|tt - 6» + lz - 7z | + m" 1 ) ' 

where we substituted 2m — 1 — I for I in the last equality. This estimate is very 
similar to the estimate of ^2 m Case 1, and hence for 7r/2 < 9 < w, 

^ 2m 

7-777 Xl* 1 < clog(TO + 1). 



2TO+ , , 



For 



Q < 2tt - 9 + 7 Z - 7; < 3tt tt 7r < 6' + 7, z -7; < 7r tt 
2 ~4 4m' 2~ 2 ~ 2 4m ' 



2m- 1 2m- 1 

1 \ - C 



and so 

2to+1 ^ 2 "2to+1 ^ (|2 7 r-6' + 7 z -7 ; |+m- 1 )(|6' + 7 z -7 J |+m- 1 ) 

2m-l 1 /)l 1 -1 

C \ - 7T — \+ TO 



2to 



+ 1 £ (k-6' + 7*+7i|+™~ 1 )(>-6'-7*-7i|+™~ 1 )' 

which is very similar to our estimate for in Case 1, and we again get a estimate 
of clogfm + 1). □ 
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5.2. Proof of Lemma [4771 

Proof. The proof of the remaining two lemmas will proceed by separating the in- 
tegral into three different regions, then dividing the sums in j and v into several 
sections, and performing estimates on each resulting section. Frequently, obtaining 
an estimate consists of bounding quotients by a constant, and then estimating sim- 
ilar "types" of sums and integrals. For the sake of brevity, we list these types here, 
and then direct the reader to the type of estimate that arises in the each piece. The 
symbols <f>\ , 4>2 and £ refer to values that are particular to the section under inves- 
tigation. For two different expressions /i and the notation {/i,/2} indicates 
that either expression satisfies that type. Finally, we note that these estimates also 
hold for sums whose range of indices are a subset of those listed below. 
Type 1: 

T^2^2 \a, Z7i d ° 



(2m + 1)2 4 {\0 - ^| , \4> 2 - 9\} + m-i ^ ^ \9 3 + £| + 
< c(log(m + l)) 2 . 



Type 2: 



Type 3: 



(2m +1)2 ^ £E(|0.+£|+ m -l)2 



-n 2 SE w, + fi +m -i d ^ 



(2m + 1)2 y 0i _ ^| , |^ - 0|} + m-i) 2 I . + CI + 

< cm log(m +1). 

Type 4: 

(2m+l)2/ 0i {|0 2 _ 0| , |0 _ ^|} + m-i £ £ (| . + £| + m -i) 2 ^ 

< cm log(m + 1). 

Type 5: 

x ,<fa i ™ 1 

(2m +1)2 7 0i {102-01,10-^H + m-i ^{|<7„(a:,j/)|, |vr - a„(.T, y)|} + r 



x 



Type 6: 
1 



<H \6j + £| +m 



(2m + l) 2 y 0i ^ {1(7,(2;, , |vr — ovfoy)!} + ™ _1 ^ \0j + CI + m" 1 
< c(log(m + l)) 2 . 
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Type 7: 

(2m +1)2 4 ^ o {|^(x,y)|, Itt-^^^D + to-i^ (|0 j+ ^| +m -i) 2 ^ 
< cmlog(m + 1). 

The above estimates are easily obtained with the following proposition and 
lemma. 

Proposition 5.1. For any real number £, 

^ 1 < j c(2m+l)log(m+l) fc=l 



^ (l^+ei+m- 1 )* " \ c(2m+l) 2 A: = 2 

Proof. By the definition of 0j, 

F r = (2m + l) fc V, ; ^— r 

m ^ 

<(2m + l) fe 2X: i - TTF , 

from which the proposition easily follows. □ 
Lemma 5.2. For (x,y) G r m , 

(5.3) > i — -, < cm log(m+l) 

^ |a v (a:, y)\ +mr l ~ 

m _^ 

(5.4) V- < cmlog(m + 1). 

Proof. Recall that (a;, y) = (rcos(0), r sin(0)). If we restrict v to < v < m/2, 
then 4> u — 4> < 7r/2, and |ov(x, y)| > |0„ — </>|, since cos(cr„(a;, y)) < cos{4> v — (f>). On 
the other hand, if <j> v — <j> > n/2, then a v (x, y) > n/2. The first inequality follows, 
since 

— 

El x ^ 1 m 2 

-, — ; n r < > t, n =- + ——< cm log(m + 1. 
v=0 Wu{x,y)\+m-^ ~ ^ +m" 1 2 tt " 

The proof of the second inequality is similar to the first. Recall that 7r — u v (x, y) = 
o- v (—x, —y), and write (—x,—y) = (r cos(<p), r sin(0)), where <f> € (7r — 4Tr ^ +2 , 7r + 
j). A similar argument shows ov(— x, —y) > tt/2 for ^ < m/2, while cr„(— x, —y) > 



4m+2 / 



|0 — </v| for m/2 < v < m, and the remainder of the proof is identical to the proof 
of the first inequality. □ 



Proposition 5.3. For < <p\ < 4>2 < n, 

r4>. 

hi 

Proof. The proposition follows from a change of variables in the integral. □ 



1 - kd e<\ cl °g(m+i) k = l 

({\e - ^\ m 2 - e\\ + m-^) k I cm k = 2 
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We introduce new notation to simplify the proof of the remaining estimates. The 
notation I„(</>i, fa) denotes the set of indices v such that fa < a v (x, y) < fa, and 
the symbol Xj(fa, fa) denotes the equivalent set of indices such that fa _! @j < fa- 

Combining A±(x,y) and A^(x,y), we obtain 



\Af(x,y)-A 1 (x,y)\ 



Are sinflj sin cr v (x, y) 



(1 - 2re ie cos(6j + cr v {x, y)) + r 2 e 2i9 )(l - 2re l9 cos(^ - a v (x, y)) + r 2 e 2i8 ) ' 

Upon substituting this into (|4.8[) and using (|5.2J) and (|5.1D , we are left with approx- 
imating 



(5.5) 



1 



EE 



i 



(2m + 1) 2 In ^ ^ sin6> + mr 1 



( sin («±& 



d 3 -(j,Ax,y) 



( sin ( g_te_feMl ) 



(18. 



In order to approximate the sine functions, the integral over [0, ir] is divided into 
integrals over three subintervals: [0, j], \j, ^1 , and 7r] . We will use the nota- 
tion Ti\, ri\, and to denote the expression (|5.5p restricted over these respective 
sub- intervals. 

Case 1: < 8 < ir/A. With this restriction on 8, the sine functions in (|5 .5[) are 
approximated by 



(5.6) 
(5.7) 

to obtain 
(5.8) H\ 



8 ±i 8j ± 2 o v (x,y) 



|0±10j±2 CT,(x,y)| 



' - C (2m+ 1) 2 7 EE + m -i 



2 

-l "i 



u=0 j=l 



1 



(|(9 + 8j + a„(x, y)\ + to" 1 ) (\9 - 8j - a v (x, y)\ + to" 1 ) 

1 

X (\8 + 8~ - a v {x, y)\ + to" 1 ) (\8 - 9~ + o v {x, y)\ + m' 1 ) ' 

First, use the inequality 8j/(8 + 8j +a„(x, y) +m~ 1 ) < 1. Considering v e X v (0, 8), 
we use the inequality 9j/(0j + 6 — cr v {x,y) + m^ 1 ) < 1, and split the sum in j into 
Tj (0, 8) and 2j (8, to obtain estimates of type 5. For (0, 7r), the inequality 

(#j + ct^(x, y) — 8 + m~ ) < 1 is used. Splitting the sum in j into Tj (0, a y (x, y)) 
and 2j (a v {x,y), ^) yields two estimates of type 3. 



Case 2: 7r/4 < 6> < 3vt/4. In this case, the factor of (sin(6>) + m -1 ) -1 in ([57 
is bounded away from zero and may be ignored. We approximate one of the sine 
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functions in (|5.5|) by 

(5.9) |sin ^ + e j+ a v ( X ,y) ^ | ^ ^ _ , _ _ ^ y)| 

and the remaining three are estimated as in (|5.6|) to obtain 

3tt 



"i r — i — 

(5.10) H\<c—^—\ VV 



fl 2 

1 



|2tt - 6» - 0j - <j„(x,y)\ + m- 1 ) (|0 - 9 3 - a v (x,y)\ + m" 1 ) 

1 



' + % - t7 w (ar,j/)| (|0- 0j + f r t/ (x,y)| 



d6>. 



First consider ^ G X„ (0, 0). Under this restriction, 2tt — — 6j — (r v (x, y) > 4p — 29 
a,nddj/(9j + — <T,j(x,y) + m~ 1 ) < 1. Splitting the sum in j in (0, #) and 2j ^) 
yields estimates of type 5. If i/ G X„ (0, 7r), first note that 0j/(9j + cr v (x,y] — 9 + 
m -1 ) < 1, and — 6j + cr v (x, y) > 2(9 — 7r/4). Substituting 2m + 1 — j for j and 
splitting the sum in j into X, 0) and lj (9, tt) yields estimates of type 5. 

Case 3: 37r/4 < 9 < ir. We may again ignore the factor of (sin(0) + m~ 1 ) 1 in 
(15. 5p . Two of the sine functions are approximated by 

(5.11) sin ( - ^ + 2 aAx > ^ ^2^-9 + 9,- *„(x, y), 

and (|5.9p , and the remaining two sine functions are approximated by (|5.6| to obtain 

„tt m m 

(5.12) H\<c - / VV( 
(2m + l) 2 73. 



m m 
4 u=0j = l 

1 



(|2tt -9-Bj- *u(x,y)\ + to" 1 ) (|0 - 0, - a v (x,y)\ + m" 1 ) 

1 

(|0 + % - a„(a;, y)| + m" 1 ) (|2tt - + 0, - a v {x, y)\ + m" 1 ) 



d0. 



First observe that 9 +a u (x, y) < 2ir, and hence 9j/(2ir — 9 + 9j — a v (x, y)+m~ 1 ) < 1. 
Considering v &X V (0, 0), use the inequality 9j/(9 + 9j — cr v (x, y) + m" 1 ) < 1, and 
substitute 2m + 1 — j for j. After splitting the sum in j into Tj (^,a u (x,y)j 
and Tj (a v (x,y),Tr), estimates of type 1 are obtained. For v G Z„ (0, tt), use the 
inequality 6^/ (0., + o~ v (x, y) — 9 + m _1 ) < 1, and split the sum in j into Ij (S, 0) 
and X, (0, 7r). The resulting expressions are type 6. □ 

5.3. Proof of Lemma 14.81 

Proof. First, we let £ = re l9 . As in the proof of Lemma (|4.7[) . we combine terms 
to obtain 

\A%(x,y) - A^(x,y)\ = 2sm9jsmcr v (x,y) |l - r 2 e 2l9 | 

x |P(C^,^(x,y))| 

(1 - re w cos(6j + a v (x, y)) + (re i0 ) 2 ) 2 (1 - re id cos(6» i - cr v (x, y)) + (re w ) 2 f ' 
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where 
(5.13) 

P(C, 9 j ,a v (x,y)) = 1 - 3C 4 - 2C 2 (3 - 2 sin 2 (9 ) - 2sin 2 (a u (x,y))) 

+ 8C 3 (1 - 2sin 2 (cr iy (a;,t/)/2) - 2 sin 2 (6» J /2) + 4 sin 2 (cr„ (z, y)/2) sin 2 (6' J /2)) 
Substituting this into (|4. 8|) and using (|5.2p and (|5.ip . it remains to estimate 



(5.14) 



(2m 



sin 



■1 />7T m m 

|P(C,^-,^(x,y)| 



( 



sin 



?+8 J +tr 1 ,(x,i/) 
2 



sin 



2 



( 



sin 



l+e ] -a v (x,y) 
2 



^j+ov(£^) 
2 



(10. 



We will again split this integral into three sub-integrals over [0, ?] , pj, ^1 , and 
[t 5 -) 71 "]) an d denote the part of (|5. 14[) associated with these subintervals by Ti.\, 
Ti. 2 , and TC^i respectively. The crucial part of the estimate is suitably approxi- 
mating \P(£,0j,(T u (x,y))\. Several different approximations will be used. These 
approximations are given in the following lemma, and are referenced as needed. 

Lemma 5.4. The function P((^,9j,a l ,(x,y)) satisfies the following inequalities. 
(El) For < 9 < tt/2, 

(5.15) \P((,9 3 ,a„(x,y))\<c (jO + Oj + m' 1 ) 2 (\9 - 6 3 \ + m" 1 ) 

+ (a„(x iy )) 2 (6 + (a v (x,y)) 2 + 9) + m- 1 )) . 
(E2) For < 9 < tt/2, 



(5.16) |P(C, ^, o-v(a:, »)| < c ((0 + m~T + (6 j + a v {x, y)Y(\9 3 - a v {x, y)\) 2 

+ (9 + m- 1 )(9 2 3 +(aAx,y)) 2 )). 
(E3) For tt/A<9 < 3tt/4, 

(5.17) |P(C, is ovfo j/)) | < c ((|0 + + m- 1 ) (|0 - 6j\ + m- 1 ) + (a v {x, y)) 2 ) . 
(E4) For tt/4 < 9 < ir, 

(5.18) \P(C,0j,<r»fav))\ 

< c ((|tt - 9 + 0j\ + m- 1 ) 2 (\tt-6- 9 3 \ + m^ 1 ) + (tt - tr^a, ?y)) 2 ) 

Proof. We first prove the estimate (I5.15[) . We define 

(5.19) Pi(C, e jt2m ) := 1 - 2C 2 (3 - 2sin 2 6 3 ) + 8C 3 cosOj - 3C 4 , 
if (C, 9 3 , a u {x, y)) := AC, 2 sin 2 a v (x, y) ± 16C 3 sin 2 



2 <iv{x,v) n 
' " 2 C0S ' 



so we may write P(£, 9j,a u (x, yj) = Pi(C, #j) + P2 (Cj 9j,(T u (x, y)). It is possible to 
factor Pi((,9j) as 

(5.20) Pi(CA) = (1 + 2CC0S6IJ -3C 2 )(1 - 2Ccos0 3 - + C 2 )- 
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The second factor of (15.20|) is approximated by (15. 2p and ()5.1|) as before, and the 
first factor may be further factored as 



( 1.21 ) (1 + 2C cos 6j - 3C) = "3(C+^(\/ 4 - sh ^ j - cos e 3 



x ( C - 3 ( V 4_sin °j +cos( 



The first factor of (|5.21|) will not be used. Using the double angle identity for 
cosine, the second factor of (|5.21[) is approximated by 



< c ( '■-'> Ml) 

< c 



4 - sin Qj + cos ( 



3 cos 9 - cos 0j — y 4 - sin 9j 



2- A/4-sin z 



Since 2 — y 4 — sin Oj < sin 9j , we obtain 

+ 6 



(5.22) Pi(C,^)<c 



Now considering P 2 (£, the double angle identities for sine and cosine are used 
to obtain 



(5.23) |P 2 -(C,^,^(x,j/)| <cK(i,j)) : 



cos 2 ' — Z COS Oj 



< c(a u (x,y)) 2 (9 + (a v (x, y)) 2 + 9 2 + m" 1 ). 



Adding the estimates (|5.22[) and (|5.23|) , we arrive at the estimate (|5.15l) . 

The proof of the estimate (|5.17p follows from replacing 9 with a constant in 
(15T5D . 

We next prove the estimate (|5.16[) . We first re- write P((, 9j,cr u (x, y)) as 



P(C, 6 s ,o v (x, y)) = 1 - 3C 4 - 6C 2 + 8C 3 + 4( 2 



sin 2 0j + sin 2 a v (x, y) 



• 2 8j -2 <7„(x,«) • 2 Oj 

.in — ^ n > — ■ — -i-sm ^ - sin 



4 ( 2 sin 2 ^ sin 



2 gv(a:,a) 
2 



+ 4(C - 1) ( 2 sin 2 & sin 2 - sin 2 U - sin 



2 gvfey) 



It is easily checked that 1 - 3( 4 - 6C 2 + 8<C 3 = -(£ - 1) 3 (3C + 1), and applying the 
double angle identity for sines, 



sin 2 9j + sin 2 a u (x, y) + 4 ( 2 sin 2 sin 



,2 cr„(x,y) 



(■20, 

[sm -f 



Finally, we may approximate the sine functions to obtain (|5.16|) . 
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Finally, we prove the estimate (15.181) . First, replace sin 2 a v (x, y) with sin 2 (7r 
<r v (x,y)) and sin 2 a " { ^ y) with 1 - sin 2 *- a ^ x ^ hi (j5T3j) to obtain 

(5.24) \P{C,e ,a v {x,y))\ < |l - 3( 4 - 4C 2 (3 - 2 sin 2 9 3 ) + 8C 3 (2sin 2 % - 1) 

+ C |^2 + (C 5 ^ I 7r - Cr ^( 2: :2/))| ■ 

The inequality \P^~((,8j,ir — <j v (x,y))\ < c(tt — a l/ (x,y)) 2 follows easily from the 
definition of P% . The first term in (|5.24p becomes |Pi(£,7r — after replacing 

3-2 sin 2 0, with 3-2 sin 2 (n-Oj) and 2 sin 2 1 with 1-2 sin 2 (^^j. Pi(C,7T-0y) 

factors as in (|5.20|) . The factor of 1 + 2 cos(9j)( — 3£ 2 can be factored further, as 

1 + 2( cos(tt - 6j) ~ 3C 2 = 1 - 2C cos 9 3 - 3( 2 

1 



-3 C 



4 - sin 



c + 



4 - sin z 



The first factor will not be used, but the absolute value of the second factor may 
be approximated by 



c + 

< c 



1 

, cos 6, 
3 * 3 



sine | 



cos ( 



< C 7T 



6 sin 



2 7T-6 



COS ( 



2 sin 2 



- 2 



< C [IT — . 



and we are able to obtain the estimate (|5.18|) . 



□ 



Case 1: < 9 < 7r/4. Approximating the sine functions in (|5.14p with (|5.6p . we 
obtain 



Hi < c 



1 



(2m + l) 2 J 



— mm 
4 \ - . 



|P(C,^,^(^2/)I 



+ 0j +<ru(x,y)\ +m^Y{\t 
1 



ovO^y)! + m X Y 



dO. 



(\9 + 9j - a u (x, y)\ + m-if (\6 - 9 3 + a v (x, y)\ + m^Y 

First considering v € X v (0, 9), note that 9 2 /(9 + 9j — <T„(x, y) + m~ 1 ) 2 < 1. Approx- 
imate l-P(C) 0j,a v (x, y)\ using (|5. 15|) . For the first term in the sum in (|5.15[) . use the 
inequality (6+0j+m ) 2 /(6+6j + a v (x,y)+m~ 1 ) 2 < 1, and then split the sum in j 
into Xj (0, 9) and X,- (6 1 , to obtain two estimates of type 7. For the second term in 
(l5~T5|) , use the inequality (9 + (cr„(x, y)) 2 + 9 2 + m^ 1 )/{9 + 9 3 + a v (x, y) + vrT 1 ) 2 < 
c(9 + m -1 ) -1 , and split the sum in j into Ij (0,9) and 2j (9, -|) to obtain two 
estimates type 4. 



For v G X v (6, it), use the inequality 8 2 /(9 3 + a v (x,y) — 9 



-1\2 



< 1, and 



approximate \P((, 0j,cr v (x, y))\ using (|5.16|) . For the first term of the sum in (|5.16p . 
use the inequality (9 + m -1 ) 2 / (9 + 9 j + cr v (x, y) + mT 1 ) 2 < 1 and then split the sum 
in j into Tj (0, <j v (x, y)) and Ij (a v (x, y), ?) to obtain two estimates of type 4. For 
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the second term of the sum in (|5.16|) . first use the fact that (9j + a^(x,y)) 2 j{9 + 
6j + a v {x, y) + to -1 ) 2 < 1, and then split the sum in j into Tj (0, a v (x, y)) and 
Tj (<t„(x, y), -|) to obtain two estimates of type 2. Finally, for the third term in 
(|5.16p , use the inequality ((a v (x, y)) 2 +9 2 )/(9 J r9 J +a lJ {x, y) +m -1 ) 2 < 1, and then 
split the sum in j into Ij (0, <r u (x, y)) and Tj (a v (x, y), ^) to obtain two estimates 
of type 4. 

Case 2: 7r/4 < 9 < 37r/4: If a v {x : y) < ir/2, the sine functions in (|5.14p may be 
approximated as in Case 1, while if a v (x,y) > 7r/2, one of the sine functions may 
be approximated by (|5.9|) and the remaining three by (|5.6|) . to obtain 




|P(CA-,M^)I 

{\9 + 9j +au{x,y)\ +m- 1 ) 2 {\9-9 ] - a v {x,y)\ +m- 1 ) 2 

1 

(\9 + ^ - o v (x, y)\ + m-i) 2 (\9 - 9 3 + a v {x, y)\ + m^ f 

m 

" {\2ir-6-6j -a„(x,y)\ + m - 1 f{\9-9 3 - a v {x,y)\ +m- 1 ) 2 

1 

X (\0 + 0j - ct u (x, y)\ + to- 1 ) 2 (\9 - 9, + <r v {x, y)\ + m^ f 

We consider first v € T v (0, n I„ (0,6*). Under this restriction, 9 + 9j + 
a u (x,y) > 7r/4 and so this factor may be ignored. Use the inequality 9 2 /(9 + 9j — 
a v (x,y) + m~ 1 ) 2 < 1, and approximate P(£,9j,a u (x,y)) by (|5.17p . After splitting 
the sum in j into Tj (0, 9) and Tj (9, %), estimates of types 2 and 7 are obtained. 

Next, we consider v € T v {9, -|) . The factor of 9 + 9j + a v (x, y) + mT 1 may again 
be ignored, and the inequality 9 2 /(9 2 + a v (x, y) — 9 + mT 1 ) 2 < 1 is used. We only 
need to estimate P((, 9j, a u (x, y)) by a constant. After splitting the sum in j into 
Ij (0, <j v (x, y)) and Tj (<j v {x, y), ^) , estimates of type 2 are obtained. 

Now we consider v £ T v (■§,#)■ For this range of v, 9 — 9j + ov(x, y) > ir/4 and 
9 2 / (9 + 9j — a v (x, y) + to -1 ) 2 < 1. We again estimate |-P(C, 9j,a u {x, y)\ by a con- 
stant. After substituting 2to+ 1 — j for j, we split the sum in j into Tj cr v {x, y)) 
and Tj {a u {x, y), n) to obtain estimates of type 2. 

To complete this case, we consider v € T v (9, tt) C\T v 7r) . Under these restric- 
tions, 9 - 9j + <r u (x,y) + mT 1 > 9 + mT 1 and 9 2 /(9 3 + a u (x,y) - 9 + mT 1 ) 2 < 1. 
Use the estimate for \P((,9j,a v (x,y)\ in (15.18|) . and substitute 2to + 1 — j for j. 
Upon splitting the sum in j into Tj and Tj (9, tt), we are left with estimates 

of types 4 and 2. 

Case 3: 37r/4 < 9 < tt. Recalling (|4.10p , a factor of tt — 9 + mT 1 is present in the 
numerator. After approximating two of the sine functions in (|5.14p by (|5.9p . and 
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(|5.11| . and the remaining two sine functions by (|5.6|) . we obtain 

-. «tt m m 

x ^CA^M 

(|2tt -6-6j- a u (x, y)\ + m" 1 ) 2 (|0 - 0j - ^(as, y)| + to" 1 ) 2 

1 

^ 

(|0 + 0j -o-„(a:,»)| +m- 1 ) 2 (|2 7 r-6' + 6i J - a v (x,y)\ + m - 1 ) 2 ' 

First, consider (0, ? ). Note that \2tt — 9 — a v (x, y) + 0j\ > tt/2, and also 

\6- cr v (x,y) +6j\ > 7r/2. Also, note that |27T - 6 - 6j - a„(x,y)\ > \tt - 9\. Ap- 
proximating \P(Ci 6j, &v{x, y))\ by a constant, estimates of type 4 are obtained. 
We next consider v £ T v ,6). We estimate \P((,6,6j,a v (x,y))\ with (|5. 18[> . 

and use the inequality 9 2 / (0 + 0j — <7„(x, y) + m -1 )" < 1. For the first term in the 
sum in (|5.18|) . use the inequality (7r+6' J -6'+m _1 ) 2 /(27r-6'-(7 y (a;, y)+6' :) +m _1 ) 2 < 1 
and split the sum in j into X, (0, tt — a v (x, yj) and X v (tt — cr v (x, y), -|) to obtain 
two estimates of type 4. For the second term of the sum in the (|5.18[) , use the 
inequality (tt — a„(x,y)) 2 /(2tt — 9 — a v (x,y) + 9j +m~ 1 ) 2 < 1 and estimates of type 
4 follow by the same splitting of the sum in j. 

Finally, we consider v £ X v (9,tt). Use the inequality 9 2 /(9j + a u (x,y) — 9 + 
mT 1 ) 2 < 1 and approximate \P(C,,Qj,cr v {x,y))\ using (|5. 18|1 . For the first sum in 
([538|) . use the inequality (tt+9 } ;-6+m~ x ) 2 /(2ir-6-a v (x, y) + 9 3 ■+rn^ 1 ) 2 < 1, and 
then split the sum in j into Ij (0, tt — 9) and Ij (tt — 9, \ ). Utilizing the factor of 
\ir— 6— Ojl+mT 1 , two estimates of type 7 are obtained. For the second term in (|5.18[) . 
use the inequality (tt— a v (x, y))(ir — 9+m~ 1 )/ (2n— 9— a u (x, y) + 9j +m~ 1 ) 2 < 1 and 
then split the sum in j into Ij (0, tt — 9) and Ij (tt — 9, ^) to obtain two estimates 
of type 7. 

□ 
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